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We use analytical methods to investigate cellular automata for traffic flow. Two differ- 
ent mean-field approaches are presented, which we call site-oriented and car-oriented, 
respectively. The car-oriented mean-field theory yields the exact fundamental diagram 
for the model with maximum velocity Vmax = 1 whereas in the site-oriented approach 
one has to take into account correlations between nearest-neighbour sites. Going beyond 
mean-field using the so-called n-cluster approach our results for Vmax > 1 are in excellent 
agreement with numerical simulations. We also present a modified cellular automaton 
which is closely related to a two-dimensional dimer model. 



1 Introduction 

In 1992 Nagel and Schreckenberg introduced a cellular automaton model describ- 
ing traffic flow of N cars on a single lane. In this model the street is divided into 
L boxes ('cells') of a certain length (for realistic applications 7.5 meters) which can 
be occupied by at most one car or be empty. The cars have an internal parameter 
('velocity') which can take on only integer values u = 0, 1, 2, . . . , Vmax- The dynam- 
ics of the model is described by update rules for the velocities and the motion of 
the cars. The update rules are given by the following four steps EJ: 

1) Acceleration: Vi v'^ = Vi + 1 (if Vi < Vmax) 

2) Slowing down: Vi — s- = — 1 (if di < Vi) 

3) Randomization: v'^ = w- — 1 (with probability p) for v[ > 

4) Car motion: Car moves v'/ sites 

Here Vi is the velocity of car i and di is the distance between cars i and i -I- 1. 
Without step 3) the dynamics would be purely deterministic and the system shows 
strong dependence on the initial condition. The randomization takes into account 
natural fluctuations in the driver's behaviour. All cars are updated simultaneously 
(parallel update) and we here use periodic boundary conditions ( "Indianapolis situ- 
ation"). For the analytical calculation it is usually advantagous to change the order 
of the steps into 2-3-4-1, since after step 1) there are no vehicles with velocity 0. 
This change in order than has to be taken into account in the calculation of the flux 
(fundamental diagram). 



pst of the results for this cellular automaton have been obtained frorp. simula- 
tions Basically there are two different methods to implement the rules B. In the 
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Figure 1: Fundamental diagram for Vmax = 5 and p = 1/2. Shown are the results of the simulations 

and the mean-field theory. 

car-oriented approach the state is characterized by the occupancy of each individual 
site (which might be empty or occupied by a car with velocity 0, 1, . . . ,Vmax)- In 
the car-oriented method on the other hand one has two lists, one with the velocities 
of the cars and one in which the distance to the next car ahead is stored. For 
small densities c ~ N/L the last method is preferable in simulations and for higher 
densities the first one. 

2 Site-Oriented Approach 

In the site-oriented approach the state of the system is characterized by the proba- 
bility Pn{t) to find at time t a site in the state n, where n = for an empty site and 
n = 1, 2, . . . , Vmax for a site occupied by a car with velocity n ? Since in mean- field 
theory one neglects correlations between n neighbouring sites one can write down 
iteration equations for the P„ in the stationary state u. These equations then allow 
a calculation of the fundamental diagram. In Fig. |l| results for the "realistic" value EJ 
Vmax = 5 are compared with the simulations. The mean- field flow is much too small 
showing the importance of correlations. 

The results of mean-fixjld theory can be improved systematically by the so-called 
n-cluster approximation era. The n-cluster method takes into account correlations 
between neighbouring sites and reduces to mean- field theory for n = 1. In order 
to get a closed set of eqiiations one uses conditional probabilities for the overlap 
of neighbouring clusters u%. Unfortunately, one then has to deal with a system of 
{vmax + 1)" nonlinear equations which in general cannot be solved analytically. 

Surprisingly, for Vmax = 1 the 2-cluster approximation already yields the exact 
result. The probabilities P{(7i, (T2) to find a 2-cluster in a state (cri, (T2) (where - due 
to the ordering 2-3-4-1 ~ <Tj = denotes an empty site and CTj = 1 a site occupied 

Remember that we do not have cars with velocity due to the ordering 2-3-4-1. 
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Figure 2: Fundamental diagram for Vmax = 1 and p = 1/2. Shown are the results of the simulations 
and the 1- and 2-cluster approximations. The 2-cluster result is exact. 



by a car with velocity 1) are given by 



P(1,0) = P(0,1) = 



1 - yi - 4gc(i - c) 

2q 



F(0,0) = l-c~P(l,0), P(l,l) = c-F(l,0), 



(1) 



where c — N/L is the density of cars and q = 1 — p. The corresponding flow is just 
/(c,p) = gF(f,0),i.e. 



f{c,p) = 



f - v/1 - 4gc(l - c) 



(2) 



Note that ~ in contrast to what one finds using random-sequential dynamics - par- 
allel dynamics yield the well-known "bunching" observed in real traffic, i.e. empty 
sites and cars attract each other (P(O)P(f ) = c(f — c) < P{1, 0)). This means that 
there exists an attraction between two cars separated by an empty site. 



For higher velocities the n-cluster approximation does not become exact reflect- 
ing the fact that long-range correlations exist for v„iax > 1 . In Fig. |^ we compare 
the results of the n-cluster approximation for n = 1, . . . , 5 with simulations. The 
n-cluster results converge quite fast (the difference between the results for n = 4 
and n = 5 is less than 1%) and already the n = 5 result is in very good agreement 
with the simulations. 



3 Car-Oriented Approach 

In the car-oriented mean-field approach one introduces probabilities Z?„ (t) for find- 
ing at time t exactly n empty sites in front of a vehicle. In this way we already have 
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Figure 3: Fundamental diagram for Vmax = 2 and p = 1/2. Shown are the results of the simulations 
and the n-cluster approximations for n = 1, 2, ... 5. 



taken into account correlations between neighbouring sites. For Vmax = 1 the sta- 
tionary solution of the equations gouverning the time evolution of these probabilities 
can be obtained quite easily u. One finds 



Do = 
D„ = 



2qc-l + ^/l - 4gc(l - c) 
p{l-Do) 



Do 

p [Do+pil-Do) 



Using this result the flux can be calculated from 

/(c,p) = (jc V D„. 



n>l 



{n > 1). 



(3) 



(4) 



In this way one recovers the result 
mean-field result is exact. 



i.e. in the car-oriented approach already the 



4 A Cellular Automaton Related to the Kasteleyn-Model 

In this Section we introduce a modified cellular automaton model in which the cars 
have no maximum velocity. Starting from the model for Vmax = 1 as described in 
Sect. 1 we modify the rules slightly. Suppose that we have applied the rules 1-4 to 
a specific car. If the car does not move in step 4 nothing is changed. But if the car 
actually drives in step 4 we carry out the car motion (by one site) and go back to 
step 2. This is repeated as long as the car actually drives in step 4. This change 
now allows cars to drive arbitrary distances (smaller than d where d is the distance 
to the next car ahead), i.e. it corresponds to a maximum velocity Vmax = W- 

It is interesting to notej-jthat this modified CA model is closely related to the 
so-called Kasteleyn model Q. This model is a two-dimensional dimer model on a 
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Figure 4: Fundamental diagram for the modified CA for different values of p. 



hexagonal lattice. By choosing the activities in an appropriate way it is possible to 
map the dimer configurations onto trajectories of a one-dimensional traffic problem 
with periodic boundary conditions EI. 

Again the fundamental diagram of this modified cellular automaton can be 
calculated exactly Already the mean-field result is exact and yields 



/(c,p) 



c(l — cjp c(l — c)v 



1 - (1 - c)p 1 
where v — p/q is the mean velocity of free traffic. 



(5) 
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